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Abstract 

Microscopic mechanism of forming the spin hquid of neutral fermions in Kondo 
lattices with nearly integer valence of the f-electrons in presented. It is shown that 
the competition between the Kondo scattering and the indirect exchange interaction 
can result in stabilization of a state with separated spin and charge excitations with 
a Fermi-type spectrum for both of them. The origin of strong antiferromagnetic 
fluctuations is briefly discussed. 
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1. The well-known Doniach dichotomy [|I| for the ground state of the Kondo lattice 
offers the antiferromagnetic state at small values of effective coupling parameter a = 
Jsf/^F and the nonmagnetic " Kondo-singlet" state at large values of a as the ground 
states of the system. Here Jgf is the on-site s/-exchange integral and ep is the energy 
characterizing the width of conduction electron band. The critical region where two 
regimes compete is determined by the condition a^Q ~ exp(— l/2a;co), or, in other terms, 
T/v ~ Tk where Tjy^K are the Neel and Kondo temperatures, respectively. We show in 
this paper that the third possibility can be reahzed at a ~ aco- the Kondo-stabilized 
spin liquid state of the RVB type is formed instead of the Kondo singlet state. This 
neutral spin-fermi liquid interacts strongly with charged conduction electrons, the low- 
temperature termodynamics is determined mainly by the spin-component of this two- 
component Fermi liquid, whereas the conduction electrons with enhanced effective masses 
are responsible for the charge transport and diamagnetic properties of the system in 
accordance with the phenomenological picture offered in ^, 

The idea of Kondo-type stabilization mechanism of the resonance valence bond state 
was offered in [^]. Coleman and Andrei described a picture of resonating valence bonds 
linked up with the electrons in a layer of Tk around the Fermi energy. However, the 
introduction of the Kondo-type mean field resulted in charge transfer from the low-energy 
conduction electrons to the neutral fermi excitations originating from the f spins. Thus 
they obtained a "charged spin liquid" like that considered in the mean- field slave-boson 
theory p. The general arguments against such a procedure were presented in and 
here we show that formation of spin liquid can prevent the system from forming the 
Kondo singlets, and the Fermi liquid with neutral spin and charged electron components 
is formed instead. 

2. We start with a standard Kondo lattice Hamiltonian H^ff which can be derived 
from the general Anderson lattice model in the case of nearly integer valence. 



Here Ek is the band level of conduction electron. Si and Si = cf^d-Cicr are the operator of 
localized f spin and conduction electron spin, respectively, a is the Pauli matrix. 

To consider this strongly coupled electron liquid in a controllable way, we start from 
the high temperature region T > Tk where the non-crossing approximation (NCA) is ac- 
ceptable and logarithmic perturbation theory is valid. We study the competition between 
the homogeneous RVB 0] and Neel ordering within the mean field approximation by 
using the Abrikosov pseudofermion operator representation (see, e.g., [§]) for the f-spins, 
Si = fia^apfip- Then the on-site and intersite contributions of the indirect RKKY inter- 



^eff - ^kC^^Cka + Jsf ^i^i " 
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action to the self energy part of the pseudofermion temperature Green's function F{ien) 
are given by the diagrams of fig. la and lb, respectively. Here the dashed and full lines 
stand for the pseudofermion and conduction electron propagators, hatched squares and 
crosshatched rectangles symbolize the sf- and ff vertex parts renormalized by the multiple 
on-site electron-spin scattering. First of these diagrams describes the standard molecular 
field due to the indirect RKKY interaction acting on the spin in a site i, and second one 
corresponds to the "anomalous" RVB-type correlations. The polarization operator in the 
diagram la is taken in a nearest neighbor approximation, 

U^itUJn) = UR{tUJn)S{q) (2) 

where 5'(q) = X^iexpiqRjj is a structure factor, Rji = \Rj — Ri\. The possibility of 
commensurate magnetic order with the antiferromagnetic vector Q such that QR = vr 
and J{R) = Jij = Jgf^Ri^) is suggested. 

We consider the renormalization of these diagrams by the on-site Kondo scattering at 
T > Tk- According to the NCA approach, this scattering processes renormalize indepen- 
dently the vertices F belonging to different cites (i,j), and this renormalization should 
be taken in the logarithmic approximation where the maximum energy variable in the 
argument of the vertex part is retained. This energy is determined by the characteristic 
frequencies of the electron Green's function G{iem) entering the polarization operator 
YlniiUn) in the RKKY interaction. 

Within the mean- field approach the transition temperatures T/v and TptvB are found 
from the self-consistent equations for the temperature Green's functions with the self- 
energy parts shown in figs.l(a,fe). These approximations are obtained in a static approx- 
imation for RKKY interaction. It was found by Doniach [jl| that the Kondo scattering 
destroys the magnetic order in a Kondo lattice at least at a > ac, and the "Kondo singlet" 
state can be formed in this strong coupling limit. Later it was recognized that one can 
expect the stabilization of RVB phase instead of the hypothetical Kondo singlet phase 
provided the inequality \Trvb — T^l -C is valid. Here we shall see that this inequality 
is realized in a critical region a- ac of the Doniach's phase diagram where Tm ^ Tk and 
the Kondo scattering nearly screens out the localized moments thus reducing the Neel 
molecular field 

B^iT) = Xj{R,T){S,) (3) 

Here Ai is the numerical factor determined by the lattice geometry. 

To find the molecular field -BAr(T) renormalized by the Kondo scattering one should 
first estimate the effective exchange integral J{T) entering the self-energy parts of Figs 
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l{a,b) (the hatched squares), 

JiR,T) =TY.G''{R,en)T''{er^,T) (4) 

n 

where the asymptotic expression for the Green's function in a coordinate representation 
for pfR ^ 1 and the spherical Fermi surface has the form 

G{R, e) = - exp(-i^ + ippR sgne) , (5) 

2'7ipfR Vp 

the vertex part 

T{e) = Jsf (l - 2a In '-f—] (6) 

\ maxje, T| J 

is the Kondo-enhanced sf-exchange integral The simple calculation give 

rep exp(——pFR) 

J{R,T)^^{ppR)a^ef de ^ (7) 

•^^ 1 + 2a In ^1 

where is the oscillatory RKKY function. 

In spite of the enhancement of the indirect exchange at T > Tk, the Kondo scattering 
eventually suppreses the AFM ordering because the same Kondo processes screen the 
localized spin entering the molecular field @. The screening factor (the insert K{T) at 
fig. lb) has the form 

ir(T) = l-2aln^/ln;^ (8) 

J-K 

within the logarithmic perturbation theory (see, e.g., ^) and formally turns into zero 
at T* = Tk{^) 1+2" > Tk- However, this function should be corrected for the exact 
solution of the on-site Kondo problem [Q, so the local spins are only partially screened 
at T approaching Tk, [K{Tk) ~ 0.37]. Nevertheless, the molecular field is diminished 
essentially in comparison with the pure Heisenberg exchange Jij. 

The spin screening is absent for the RVB molecular field Brvb- To find this field one 
should calculate the self-energy part of Fig. lb in a static approximation Ilq{iuJn) nq(0). 
Then the field Bfcvb is given by the equation (see, e.g., |^, [T0| ) 



BRyB{T) = \2J{R,T)^RVB (9) 

where ^rvb = J2a < /it /jo- is the "anomalous" RVB correlator. The homoge- 
neous RVB state possesses the fermi-type excitation spectrum with a dispersion u{k) = 



4 



A2 J(-R, T)ARVBS(k) To find the mean-field transition temperature Trvb one should solve 
a self-consistent equation 



^RVB = [zN) 2^ S (k) tanh — (10) 

k ^-^ 

(cf. 0). Here z is the coordination number. If the lattice geometry is such that several 
coordination spheres (m) give negative sign and comparable values, then 

the RVB molecular field has the form |T^ 

BRVBiT) = \,AnvBT.ARm,T)S^"'\k) . (11) 

m 

where m is the number of coordination sphere. To simplify the numerical calculations 
we used the approximation Brvb{T) ~ \2J{Ri,T)AjiVBS^^\\i) where A2 is a numerical 
factor. 

Then, the molecular fields Bn(T) and Bf>vB{T) are given by 

Bn{T) ^ B%{T)K{T/Tk) ln-"(T/Tx), (12) 
BnvBiT)c:^B%y^iT)\n~-iT/TK). 

where index (0) stands for stands for a "pure" Heisenberg interaction Jij without Kondo 
renormalization. In the convolution procedure of eqs.(^,(|^) the logarithmic enhancement 
is partially "integrated out", and n is close to 1 for the reasonable values of the model 
parameters. When comparing these equations we find that exchange enhancement favors 
both Neel and RVB states but the spin screening infiuences only Neel field thus weakening 
the trend to the AFM ordering. The possible contribution of several coordination spheres 
in Brvb (0) is an additional factor stabilizing the spin-liquid state, because the staggered 
Neel field possesses the definite antiferromagnetic vector Q. We have taken this possibility 
into account by varying the factor A2 in Brvb- 

The phase diagram (T^,Tjiyb vs a) depends also on the form of the electron Fermi 
surface which infiuences the coordinate dependence of the function ^{x). We calculated 
the transition temperatures for the spherical and cylindrical Fermi surfaces with $(x) ~ 
nx~^ cos2x and ^{x) ~ —2x~'^ sin 2a:, respectively, and found that the region of existence 
of the RVB phase is larger in the latter case which is close to the real situation in CeRu2Si2 

0- 

The example of numerical solution of the mean field equation for TN,TFtvB with cylin- 
drical Fermi surface is presented in figure 2. Here again the index "0" marks the transition 
temperatures for the bare Heisenberg interaction, T^,T]jy^ ~ epa^^iPFR) 
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It is known that > T'^y^ in a 3D case and the RVB state does not arise in three 
dimensions, but here we find the the Kondo-processes change essentially the Doniach's 
diagram (T^, Tk vs a) in a critical region a ~ etc- We have instead a modified diagram 
(Tiy,T^yB vs «) with a wide region of spin liquid phase overcoming the magnetic phase 
and arising instead of the Kondo-singlet phase. 

The critical region is characterized by the temperature sequence Tjivb > > Tk, 
all of them being of the same order. This means, firstly, that the spins entering the RVB 
pairs are screened essentially by the Kondo interaction, and, secondly, that the neutral 
spin liquid is close to the antiferromagnetic instability. On the other hand, the Kondo 
scattering is "quenched" at T T^vb > Tk- This means that the Kondo singlet state is 
not formed in this case, so the Kondo mean field (cffi) which could glue together the spin 
and charge excitations |]^, Q equals zero. Then, the Kondo temperature is no more the 
singular point of the theory, and we can transit to the low-temperature region T < Tk 
where the two-component Fermi liquid occurs. This liquid contains the slow electrons 
with the energy e < Tk which interact strongly with the neutral spin fermions with the 
energy scaled by T^yb- The coupling constant can be estimated as ~ J{Trvb)- Just this 
picture was offered in 0, ^ ^ for description of the heavy fermion state in the rare-earth 
intermetallides with nearly integer valence (see also |jl2|, p!3| . 

3. In conclusion, we obtained that the Kondo processes in Kondo lattices at T above 
Tk prevent the system from forming the Neel state in a critical region a ~ of the 
Doniach diagram. The same processes stabilize the spin liquid state which, in turn, 
results in binding the spins in the RVB pairs and quenching the Kondo processes at 
T ~ Trvb- At low temperature we have a two-component Fermi liquid of strongly 
interacting spin and charge fermions instead of the Kondo singlet, and Tryb plays the 
role of characteristic temperature of the heavy fermion behavior. However, since Tbvb is 
close to T/v, the AFM fiuctuations should be strong in this spin liquid, and the magnetic 
order can be restored under small external infiuences. 

There are many essential effects which are left beyond the mean-field approximation. 
First of all, the mean-field approximation introducing the average ^rvb violates the 
gauge invariance of the Hamiltonian (|I]) (see, e.g., IQ). The fiuctuations of the phase of 
pseudofermions should result in transformation of an abrupt phase transition described 
by equation (|TU|) into a smooth crossover. It is known [0, |^ that the homogeneous 
fiuctuations of the gauge field in a 2D t — J model (chirality fiuctuations) result in infrared 
singularity of the gauge field fluctuation correlator. In a 3D sf-exchange model the RVB 
transition is close to the antiferromagnetic instability according to the phase diagram of 
fig. 2, hence the phase fiuctuations should be essentially inhomogeneous. Thus, we believe 
that the infiuence of the gauge fiuctuations is not so destructive for the RVB pairs as in 
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the 2D case, and the spin-charge separation persists in the heavy fermion systems at low 
temperatures. 

When the retardation effects (the frequency dependence of the polarization operator) 
are taken into account in the diagram lb, the "polaron renormalization" ~ oP'uj In ^ and 
the damping ~ o?T of the spin-fermion spectrum arise at T < T^y^- The spin fermions, 
in turn, should influence noticeably the electron spectrum p|, |^, [1^, 0. One should hope. 



however, that these processes will not change in principle the picture of forming the two- 
component quantum liquid consisting of charged electrons and neutral spin fermions 
with different effective masses which share the responsibility for the thermodynamic and 
transport properties of the heavy fermion systems. 
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Figure captions 

Fig.l. Self-energy parts for the pseudofermion Green's functions Ga (a) and (b). 

Fig. 2. Modified Doniach diagram for competing Neel and RVB phases calculated with 
the values of 2; = 6, A1/A2 = 2.1, p^R — 2.88. ctco and etc are the critical points where 
the AFM solution disappears in Doniach's and modified Doniach's pictures, respectively. 
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